Diversity of discrete breathers observed in a Josephson ladder 
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We generate and observe discrete rotobreathers in Josephson junction ladders with open boundaries. 
Rotobreathers are locahzed excitations that persist under the action of a spatially uniform force. 
We find a rich variety of stable dynamic states including pure symmetric, pure asymmetric, and 
mixed states. The parameter range where the discrete breathers are observed in our experiment is 
Umited by retrapping due to dissipation. 
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Nonlinearity and lattice discreteness lead to a generic 
class of excitations that are spatially localized on the 
scale comparable to the lattice constant. These excita- 
tions, also known as discrete breathers^ have recently- 
tracted a lot of interest in theory of nonlinear lattices 
It is believed that discrete breathers might play an impor- 
tant role in the dynamics of various physical systems con- 
sisting of coupled nonlinear oscillators. It has been even 
said that that for discrete nonlinear systems breathers 
might be as important as solitons are for the continuous 
nonlinear media. 

There have been several recent experiments that 
reported on generation and detection of discrete 
breathers in diverse systems. These arja low-dimensional 
crystala3, anti-ferromagnetic materialsQ, counkd optical 
waveguideaj, and Josephson junction arrayajfl. By us- 
ing the method of low temperature scanning laser mi- 
croscopy, we have recently reported the first direct visu- 
alization of discrete breathersQ. 

In this paper we present new measurements of localized 
modes in Josephson ladders. Using the same method as 
in our first experiment, we study an even more tightly 
coupled lattice of Josephson junctions and observe a rich 
diversity of localized excitations that persist under the 
action of a spatially uniform force. 

A biased Josephson junction behaves very similar to its 
mechanical analog that is a forced and damped pendu- 
lum. An electric bias current flowing across the junction 
is analogous to a torque applied to the pendulum. The 
maximum torque that the pendulum can sustain and re- 
main static corresponds to the critical current 1^ of the 
junction. For low damping and bias below /c, the junc- 
tion allows for two states: the superconducting (static) 
state and the resistive (rotating) state. The phase differ- 
ence ip of the macroscopic wave functions of the super- 
conducting islands on both sides of the junction plays the 
role of an angle coordinate of the pendulum. According 
to the Josephson relation, a junction in a rotating state 

generates dc voltage V = where (...) is the 

time average. By connecting many Josephson junctions 
by superconducting leads one gets an array of coupled 
nonlinear oscillators. 

We perform experiments with a particular type of 



Josephson junpjtkun array called Josephson ladder. Theo- 
retical studiesErty of these systems have predicted the 
existence of spatially localized excitations called roto- 
breathers. Rotobreathers are 27r-periodic solutions in 
time that are exponentially localized in space. 
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FIG. 1. Optical (a) and schematic (b) view of a linear lad- 
der, (c) Spatially homogeneous whirling state measured by 
the laser scanning technique. 

Measurements are performed on linear ladders (with 
open boundaries) consisting of Nb/Al-AlOa^/Nb under- 
damped Josephson tunnel junctionstll. An optical im- 
age Fig. |^(a) shows a linear ladder that is schematically 
sketched in Fig. Q(b). Each cell contains 4 small Joseph- 
son junctions. The size of the hole between the supercon- 
ducting electrodes which form the cell is about 3x3 jim^ . 
The distance between the Josephson junctions is about 24 
/im. The number of cells N in the ladder is 10. Measure- 
ments presented in this paper have been performed at 5.2 
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K. The bias current Ib was uniformly injected at every 
node via thin-film resistors Rb = 32 Q. Here we define 
vertical junctions (JJy) as those in the direction of the 
external bias current, and horizontal junctions (JJh) as 
those transverse to the bias. The ladder voltage was read 
across the central vertical junction. The damping coeffi- 
cient a = ^^o/{27rIcCR'^g) is the same for all junctions 
as their capacitance C and sub-gap resistance R^g scale 
with the area and Cn/Cy = RsgV / RsgR- The damping a 
in the experiment can be controlled by temperature and 
its typical values vary between 0.1 and 0.02. 

There are two types of coupling between cells in a 
ladder. The first is the inductive coupling between the 
cells that is expressed by the self-inductance parameter 
Pi, = 27rL/cv/^o, where L is the self-inductance of the 
elementary cell. The second is the nonlinear Josephson 
coupling via horizontal junctions. The ratio of the hori- 
zontal and vertical junction areas is called the anisotropy 
factor and can be expressed in terms of the junction crit- 
ical currents 77 = /ch/^cV- If this factor is equal to zero, 
the vertical junctions are decoupled and operate inde- 
pendently one from another. On the other hand, if this 
factor goes to infinity the ladder behaves like a paral- 
lel one-dimensional array and no rotobreather can exist 
since no magnetic fiux can enter through the horizon- 
tal junctions. Thus, it is an important challenge to in- 
crease the anisotropy as far as possible. In this work we 
present measurements for the highest anisotropy factor 
(77 = 0.56) studied up to now— In contradiction to the ex- 
isting theoretical predict ionE^I for this anisotropy value, 
in the studied parameter range with moderate dissipation 
we find a rich diversity of localized excitations. 

To briefiy introduce the role of parameters we like to 
quote the equations of motion for our system (see Ref. |^ 
for details): 

V>7 + a^Y + sm^Y = l- ^i-^V>7 + V<^f_i - ^V>f-i), 

where (pY , (/pp, (^p are the phase differences across the Ith 
vertical junction and its right upper and lower horizontal 
neighbors, \/(pi = (pi^i - (pi, A(pi = (pi^i + (pi-i - 2(pi, 
and 7 = Ib/IcV is the normalized bias current. 

In order to generate a discrete breather in a ladder we 
used the technique described in Ref. |^. We used two extra 
bias leads for the middle vertical Josephson junction (cf. 
Fig. |l|(b)) to apply a local current /locai > ^v- This cur- 
rent switches the vertical junction in the resistive state. 
At the same time, being forced by magnetic ffux conser- 
vation, one or two horizontal junctions on both sides of 
the vertical junctions also switch to the resistive state. 
After that /locai is reduced and, simultaneously, the uni- 
form bias Ib is tuned up. In the final state, we keep 
the bias Ib smaller than /cv and have reduced /locai to 
zero. By changing the starting value of /locai it is possible 



to get more than one vertical junction rotating. When 
trying the slightly different current sweep sequence de- 
scribed in Ref. we generated mainly states with many 
rotating vertical junctions. 

In all measurements presented below, we have mea- 
sured the dc voltage across the middle vertical junc- 
tion as a function of the homogeneous bias current Ib- 
We used the, method of low temperature scanning laser 
microscopytd to obtain electrical images from the dy- 
namic states of the ladder. The laser beam locally heats 
the sample and changes the dissipation in an area of 
several micrometers in diameter. If the junction at the 
heated spot is in the resistive state, a voltage change will 
be measured. By scanning the laser beam over the whole 
ladder we can visualize the rotating junctions. 
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FIG. 2. Various localized states (discrete rotobreathers) 
measured by the low temperature scanning laser microscope: 
(a)-(d) asymmetric rotobreathers, (e)-(li) symmetric roto- 
breathers. Region M is illustrated in Fig. hJ(b). 
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Various measured ladder states are shown in Fig. g 
as 2D gray scale maps. The gray scale corresponds to 
the measured voltage response during the laser scanning. 
In Fig. ^(a) we present the simplest of observed states, 
an asymmetric single-site rotobreather, where one verti- 
cal Josephson junction and the upper adjacent horizon- 
tal Josephson junctions are in the resistive state. On 
the right side of the plot we show the corresponding 
schematic view with arrows marking the rotating junc- 
tions. In the case of Fig. ||(b) two vertical junctions and 
two horizontal junctions are rotating that corresponds to 
an asymmetric two-site breather. An asymmetric three- 
site breather and an asymmetric four-site breather are 
shown in Fig. |^(c) and (d) respectively. Here, the voltage 
of the whirling horizontal and vertical junctions are equal 
due to magnetic flux conservation. The magnetic flux 
enters the ladder through one horizontal junction, goes 
through the vertical ones and leaves the ladder through 
another horizontal junction. When a single magnetic 
flux quantum is passing through a Josephson junction 
its phase is changing by 27r. In contrjast to our previ- 
ous measurements of an annular ladderLl,here in ladders 
with open boundaries we observe asymmetric breathers 
as frequently as symmetric breathers. The Ib-V charac- 
teristics of these localized states is presented in Fig. ^(a). 
Particular states indicated on the plot are stable along 
the measured curves. The more junctions are whirling, 
the higher is the measured resistance. 

Another type of discrete breathers observed in our ex- 
periment is shown in the Fig. ^(e)-(h). Figure ||(e) il- 
lustrates a symmetric single-site breather with one ver- 
tical and all four adjacent horizontal junctions are in 
the resistive state. Though we call this states symmet- 
ric, the upper and the lower horizontal junctions may 
in general have different voltages. The simplest voltage- 
symmetric case would be when each horizontal junction 
voltage is equal to half of the voltage of the whirling ver- 
tical junctionQ. We also observed a variety of multi-site 
breathers of this type. A two-site, a three-site, and a 
four-site breather is shown in Fig. |2| (f), (g), and (h), re- 
spectively. The corresponding Ib-V characteristics over 
the stability range of symmetric breathers are presented 
in Fig. 1(b). 

As mentioned above. Figs, ^(a) and (b) show the Ib-V 
characteristics with /locai = of asymmetric and sym- 
metric rotobreathers, respectively. The voltage V is al- 
ways recorded locally on the middle vertical junction 
which was initially excited by the local current injec- 
tion. The vertical line on the left side corresponds to the 
superconducting (static) state. The rightmost (also the 
bottom) curve accounts for the spatially uniform whirling 
state (all vertical junctions rotate synchronously and hor- 
izontal junctions are not rotating). An electrical image 
of this state is shown in Fig. |l|(c). The upper branches in 
Fig. H represent various localized states. The uppermost 
branch corresponds to an single-site breather, the next 
lower branch to a two-site breather, and so on. 
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FIG. 3. Ib-V characteristics for asymmetric (a) and sym- 
metric (b) breathers in a ladder with the parameters N = 10, 
IcV = 320 /iA, T] = 0.56, and /3l = 4.3. The grey region indi- 
cates the frequency range of the upper plasmon band. Up is 
the plasma frequency. 

It is easy to show that the voltage V^ at which the verti- 
cal and horizontal junctions switch back to the supercon- 
ducting state is the same. By using Rsgv/RsgH = V 
Icu/IcV = ^ we get with 1^ oc Ic that V^ = RsgvIrV = 
^sgH^rH- For the uniform state we observed V^ ~ 0.6 
mV. At about the same voltage, the whirling horizontal 
and vertical junctions for the asymmetric breather are 
retrapped to the static state as can be seen in Fig. ||(a). 
The symmetric breathers show different behavior. As- 
suming Vu ^ Vv/2, we can expect that already at the 
voltage V = 2Vy ^ 1.2 mV the horizontal junctions 
should trap into the superconducting state. But if the 
voltages of the top and bottom horizontal junctions are 
not equal while decreasing /b, the retrapping current in 
the junction with the lower voltage is reached earlier and 
the retrapping occurs at a higher measured voltage. 
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FIG. 4. The retrapping current measured (points) and cal- 
culated (line) by Eq. ^ for the asymmetric (a) (V = K) and 
the symmetric (b) {V = 2Vr) rotobreathers. The insets show 
the reduced electric circuits for each case. 

We found that in order to explain the retrapping cur- 
rent of the breather states the bias resistors Rb have to 
be taken into account. Initially, these resistors have been 
designed for providing a uniform bias current distribution 
in the ladder. Assuming that the Josephson junctions in 
the superconducting state are just short circuits we get 
the electric circuits which are shown as insets in Fig. |^. 
The total current injected in the breather region is cal- 
culated by using the resistance R of vertical junctions 
measured in the spatially-uniform state. Thus we derive 
the retrapping current as 



(iV+l)/B 



(A/- - M + 1) (AT + 1)[M + (2 - S)r]] 
(1 + S)Rb ^ 



MR 



(4) 



where M is the number of whirling vertical junctions. 
The parameter 6 is equal to for asymmetric breathers 
and to 1 for symmetric breathers. The calculated retrap- 



ping current is compared with the experiment in Fig. |^. 
For the asymmetric breathers the agreement is fairly 
good. We believe that the larger discrepancy for the 
symmetric breathers is due to the assumption Vu = Vy /2 
that we imposed in derivation of the retrapping current. 
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FIG. 5. Ib-V characteristics of three-site breathers. 

One interesting feature in our experiment is the dy- 
namical behavior of the local states in the linear ladder 
which differs-, from the previous observations in the an- 
nular ladderQ with the anisotropy parameter 77 = 0.44. 
In Fig. H this feature is illustrated for 3-site breathers. 
After the creation of a symmetric rotobreather (middle 
curve) we are lowering the uniform bias current /b until 
the retrapping current of one of the horizontal junctions 
at each side is reached. Here we observe switching to 
the corresponding asymmetric 3-site rotobreather state. 
By further lowering of the bias current we get to the 
point where both the vertical and the horizontal junc- 
tions reach their retrapping currents and the whole lad- 
der goes into the superconducting state. If, instead of 
lowering the bias current, we start increasing it another 
switching point is observed where the previous symmet- 
ric breather state is recovgred. In contrast to this behav- 
ior, in the annular ladderEl the observed lower instability 
of the symmetric breather led to an increase of voltage. 
Moreover, this switching turned the number of rotating 
junctions to increase, with few more vertical junctions 
starting to rotate. 

In addition to the simplest hierarchic states described 
above we also found a large variety of more complex lo- 
calized dynamic states. Some examples are presented in 
Fig. ^. Fig. 11(a)- (c) are sections from the middle re- 
gion (M) [cf. Fig. [^(b)] of the measured linear ladder. 
Fig. ^(a) shows an asymmetric 6-site breather with top 
and bottom horizontal junctions whirling on its sides. 
Fig. ^(b) illustrates a 3-site rotobreather which is sym- 
metric on one side and asymmetric on the other. Another 
very peculiar localized state is presented in Fig. 0(c). 
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The rightmost vertical junction is rotating with a lower 
frequency then the other four vertical junctions on the 
left. This can be understood from the appearance of 
an extra rotating horizontal junction in the interior of 
this state. We can interpret it as a 4-site breather cou- 
pled to a single-site breather. All these localized states 
would be topologically forbidden in the case of an an- 
nular ladder as the magnetic flux inside the supercon- 
ducting circuit should not accumulate, but they are not 
forbidden in a ladder with open boundaries. The last two 
pictures (d) and (e) are taken near the border of the lad- 
der. We find here truncated asymmetric and symmetric 
6-site breathers. The marginal vertical junction does not 
require any horizontal junctions to rotate. 
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FIG. 6. More complex non-uniform states measured in the 
middle (M) and on the border (B) regions of the ladder. 
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In summary, we presented observation of a large vari- 
ety of spatially-localized dynamic rotobreather states in 
Josephson ladders with open boundaries. These states 
can be excited by the local current injection and sup- 
ported by the uniform current bias. We observe both 
symmetric and asymmetric rotobreather states predicted 
m Ref. [0[ as well as more complex mixed states. We 
believe that the region of existence and stability of ro- 
tobreathers does sensitively depend on dissipation, dis- 
creteness, and anisotropy parameter of the ladder. 
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